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SPACE-BASED GRAVITATIONAL WAVES AND LISA BAYESIAN ANALYSIS FOR GW 
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e LISA (Laser Interferometer Space Anten- 
nae) is a near-future mission ( 2034) to 
measure gravitational waves (GW). 

e Emitted from some of the most elusive and 


GALACTIC BINARIES FOR CALIBRATION 


LISA is expected to detect around 2.5 x 104 
WDB systems in the galaxy, around 107 of 
which will have EM observations from Gaia and 
the LSST [3]. Their gravitational waveforms 
are nearly continuous, and therefore considered 
verification binaries since their parameters can 
be compared to those obtained from EM analy- 
ses. This means they provide a calibration sys- 
tem independent of those implemented in mis- 
sion design. 


h(t) = Asin[2n(f + £)t + @p(t) + go] 
Allows for: 


e Estimation of temporal gaps expected to 
be present in data by simple phase shift 
estimation in waveforms (tt expected to 


be small, but further studies may examine 
this). 


e Estimation of total error in 
phase (69) and amplitude (6A) 
as additional hyper-parameters in 
MCMC. 


energetic events in the Universe. 

e Accessible source frequency range 10~° to 
10°! Hz (unlike ground-based instruments 
such as LIGO) 

e Measurement test masses aboard three 
separate spacecraft in heliocentric orbit. 

e Separated by a distance of L = 2.5 x 10° km. 

e Each spacecraft with two optical links for 
digital interferometry. 

e GW detected by change in test mass sep- 
aration 6L like LIGO. However many dif- 
ferences and challenges in the data anal- 
ysis methods for spaced-based GW detec- 
tors are present and will be discussed here. 


TIME-DELAY RANGING (TDIR) INMCMC 


Unlike ground-based LIGO, the distances be- 
tween test masses are not fixed, leading the vari- 
able arm-length of the interferometers to cause 
laser frequency noise (LFN) to be present in the 
data; amounting to 10° times greater strain (34) 
than the expected strain due to gravitational 
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=). Time-delay interferome- 


try (TDI) is the accepted method of suppress- 
ing LEN by delaying the phase measurements 
transmitted from individual spacecraft by a du- 
ration D corresponding to the relative separa- 
tion, then combining delays to cancel the LFN 
which requires 100 ns timing accuracy. 
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Rather than delaying the signals in pre- 
processing,  fractional-delay interpolation 
(FDI) re-samples the signal with the non- 
integer delay D, given equally-spaced raw 
phase measurements. The delay will then be 
estimated as an additional parameter in the 
same MCMC used in GW signal analysis. 
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A currently implemented analysis algorithm 
Bayeswave [4] used in ground-based GW infer- 
ence will form the foundation for adaptation in 
space-based (SB) measurements. 


p(hls, M) = Se ae 

e Infers posterior distribution of waveform 
h given data and model, p(h|s, M) 

e Model is either: GW signal, noise, or 
glitch. 

e Multi-dimensional reversejump Monte 
Carlo Markov Chain (MCMC) in 
Bayeswave finds evidence and complexity 
(number of wavelets for reconstruction) 
for each model to determine confidence 
in astrophysical origin and estimates 
signal waveform parameters for source 
characterization. 
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This method has been consistent in all of the 
last three LIGO BBH detections. A similar al- 
gorithm for SB detection with LISA will be de- 
veloped and tested on LISA Pathfinder data to 
determine how noise models may need modifi- 
cation due to different noise sources and prop- 
erties. 


NEXT STEPS FOR FDI 


FDI convolves the data with a sine cardinal fi- 
nite impulse response (FIR) with a defined de- 
lay [2]. 


s(n—D) = s(n) * sinc(n— D) 


Since TDI requires 100 ns timing accuracy, the 
FIR approximation yielding the least interpola- 
tion error without loss of data will be found. 


~2infD 
e=max(|H(f)-e *  |impz<f<1Hz) 


Sinc filters currently being tested include Black- 
mann and Hammond windowed-sinc filters 
and the LaGrange filter. Of these, the Black- 
mann windowed-sinc filter w? is expected to be 
optimal in testing the relative difference in the 
delayed re-sampling of the data with original 
simulated data. 


Wp = 0.42 + 0.5 cos(q=7) + 0.08 cos( 7) 
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